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Long distance electron transfer (ET) plays an important part in many biological processes. Also, 
fundamental understanding of ET processes could give grounds for designing miniaturized electronic 
devices. So far experimental data on the ET mostly concern ET rates which characterizes ET 
processes in whole. Here, we develop a different approach which could provide more information 
about intrinsic characteristics of the long-range intramolecular ET. A starting point of the studies 
is an obvious resemblance between ET processes and electric transport through molecular wires 
placed between metallic contacts. Accordingly, the theory of electronic transport through molecular 
wires is applied to analyze characteristics of a long-range electron transfer through molecular bridges. 
Assuming a coherent electron tunneling to be a predominant mechanism of ET at low temperatures, 
it is shown that low-temperature current-voltage characteristics could exhibit a special structure, 
and the latter contains information concerning intrinsic features of the intramolecular ET. Using the 
Buttiker dephasing model within the scattering matrix formalism we analyze the effect of dephasing 
on the electron transmission function and current-voltage curves. 

PACS numbers: 05.45.+b, 02.60.+y 



I. INTRODUCTION AND BACKGROUND 

It is well known that long distance electron transfer 
(ET) plays an important part in many biological pro- 
cesses such as signal transduction across membranes, 
photosynthesis, enzyme cathalysis and other reactions 
in biological cells [0, Q- A fundamental understand- 
ing of these reactions also has substantial implications 
for further miniaturization of electronic devices (3, 0, ■ 
Theoretical and experimental investigation of ET pro- 
cesses in biological macromolecules such as proteins and 
DNA lasts more than three decades but it still remains a 
very active area of research. ET kinetic data on macro- 
molecules are increasingly available (for recent experi- 
mental results, see e.g. ja0 an d references therein). 

It has been established that molecular ET is essen- 
tially a combination of nuclear environment fluctuations 
and electron tunneling. Due to the large distances be- 
tween donor and acceptor subsystems in macromolecules 
(~ lnm), a direct electronic coupling between them is 
extremely small, so the ET is provided by intervening 
molecular bridges, giving rise to sets of intermediate 
states for the electron tunneling. The expression for the 
ET rate including both electronic and nuclear factors was 
first proposed by Marcus 0, 0, |HJ and can be written as 
follows: 

K ET = K el K n v n . (1) 

Here, K e i is the electron transmission coefficient, K n is 
the nuclear transition or Franck-Condon factor, and v n 
is the effective nuclear vibration frequency. 

The enormous size of biological molecules makes 
straightforward calculations of K e i extremely difficult. 
As a result of numerous theoretical studies, some meth- 
ods have been developed to simplify the analysis of ET 



processes. These methods are mostly based on a general 
model of prim ar y or hig hly prob able pathways (see e.g. 

Refs. mmElQEHEllIl)- The P° int of this ap- 
proach is that an electron involved in a long range ET 

process, with a high probability follows a few primary 
tunneling pathways moving through the macromoleculc. 
Therefore large portions of the original macromoleculc 
do not significantly contribute to the electron transfer, 
and can be excluded from further calculations. Keep- 
ing only these primary pathways, one can reduce original 
molecular bridges to much simpler chain-like structures. 
This enables to develop a numerical analysis of elecron 
transport along the chains. 

Assuming the coherent tunneling (superexchangc 
mechanism) to predominate one can treat a long-range 
ET process as a sequence of tunnelings between orbitals 
of the considered donor-bridge-acceptor system. The lat- 
ter are sites of the system, and any atom included into 
consideration is represented as a set of sites correspond- 
ing to its states associated with valence electrons. The 
electron transmission coefficient k e i equals |Tda| 2 , where 
Tda is the transition amplitude between donor and ac- 
ceptor subsystems. The transition amplitude can be pre- 
sented in the form: 

j,mn _ j,0^ _|_ ^2 V>,„, C',,., I, „. (2) 

ij 

Here, T^„ is the direct coupling of the donor site "m" 
to the acceptor site "n" which is usually weak in macro- 
molecules and therefore can be neglected. The summa- 
tion in the second term is carried out over the bridge 
sites involved in the ET process; T> m i and Aj n are the 
couplings of the donor/acceptor to the bridge; G,j are 
matrix elements of the Green's function between bridge 
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sites i and j: 

G l3 =<i\{E-H)- l \ ] > (3) 

where, H is the Hamiltonian of the bridge, and E is the 
tunneling energy. 

This method of calculation of Tda was repeatedly em- 
ployed in numerous papers |l8j . Obtained results were 
succesfully used to justify semiempirical expressions for 
k e i which were applied for interpretation of experimental 
results, such as the pathway model (see Ref. |(|), or the 
average packing density model [llij . 

Thus far, insufficient attention was paid to systemati- 
cally analyze the contribution of donor/acceptor coupling 
to the bridge to the long-range ET processes. Usually, in 
theoretical studies concerning ET processes both donor 
and acceptor are presented as single sites, whereas in re- 
alistic biological macromolecules the donor/acceptor sub- 
system often has a complex structure and includes a set 
of sites providing effective coupling to the bridge. Cor- 
respondingly, the bridge has a set of entrances and a set 
of exits which an electron involved in the ET process can 
employ. At different values of the tunnel energy E dif- 
ferent sites of donor and acceptor subsystems can give 
predominant contributions to the ET. In other words, at 
different values of E preferred entrances to the bridge and 
exits from it, as well as primary pathways for electrons, 
differ. 

ET processes in complex molecular systems occur 
in a highly disordered nuclear/solvation enviromcnt. 
Stochastic fluctuations of the environment electric po- 
tential cause the electronic phase-breaking effect. Also, 
incoherent component in the electron transport could 
originate from the electron scattering on intramolecu- 
lar vibrational modes |2^| and due to some other dis- 
sipative processes Q- The dephasing could strongly 
affect the character of the ET processes destroying elec- 
tron pathways in the molecules and providing a transition 
from coherent superexchange to the completely incoher- 
ent sequ ental hopping mechanism of the electron transfer 
jH EE El El The effect of inelastic scattering in the 
ET processes is often described within the Buttiker ap- 
proach . The electron propagation along the molecu- 
lar bridge is considered as a multichannel scattering prob- 
lem which is treated u sing the scattering matrix formal- 
ism El|^|j^|j^|33|,|jM|]. Also, some other formalisms 
are employed to treat the inelastic scattering effects on 
the ET characteristics which are discussed elsewhere ■ 

ET processes reveal a strong resemblance with the elec- 
tron conductance through molecular wires |36| . This 
resemblance is a starting point for the present work. 
We believe that studies of electronic transport through 
macromolecules can provide important information con- 
cerning quantum dynamics of electrons participating in 
the ET processes. In this paper a simple approach is 
proposed which shows that some intrinsic characteristics 



of the intramolecular ET such as pathways of tunnel- 
ing electrons and distinctive features of donor/acceptor 
coupling to the bridge at different values of tunnel en- 
ergy can be obtained in experiments on low-temperature 
electrical transport through the corresponding molecules. 
These data may be available under moderate electronic 
dephasing |2jfl , when the structure of electron transmis- 
sion function containing the desired information is not 
completely washed out. A similar approach could also 
be developed to analyze characteristics of intermolecular 
ET processes in extended bio-systems. 

The organization of this paper is as follows. In Sec. 
II the adopted model is introduced to analyze coherent 
electron tunneling through a molecular bridge. Within 
this model we take into account complex structure of 
donor and acceptor subsystems which is typical for re- 
alistic macromolecules. In the Sec. Ill we analyze how 
electronic scattering influences the transmission function. 
We follow Buttiker dephasing model, and we apply it to a 
simple single-site bridge. This enables to solve the prob- 
lem mostly analytically. In the last Section we discuss 
how these results could be used to get additional informa- 
tion concerning intrinsic characteristics of ET processes. 

II. THE COHERENT ELECTRON TRANSPORT 

Existing theories of electrical transport through macro- 
molecules are mostly based on a very simple model 
simulating both electrodes by semi-infinite tight-binding 
chains attached to the ends of the molecule. The latter 
is also simulated as a single tight-binding chain of sites 
(see e.g. [2Jl|3^|3lJ,|3^|3^|33,|39|). Here we adopt a re- 
sembling model to analyze the electric transport through 
realistic macromolecules including donor, bridge, and ac- 
ceptor subsystems, but we take into account that donor 
and acceptor subsystems are complex and include many 
sites which exhibit a significant coupling to the bridge. 

To simplify first steps of further calculations we sep- 
arate out the effect of electron dephasing, assuming co- 
herent electron tunneling through the macromolecule to 
be the predominant mechanism of electrical transport. 
Correspondingly, we treat the electron transport through 
the molecule as a sequence of tunnelings between poten- 
tial wells. Each well represents one atomic orbital of the 
donor-bridge-acceptor system. An electron which partic- 
ipates in the transport is treated as moving in an array 
of these potential wells which are the sites of the array. 
Within this approach, any atom is represented by a set 
of sites corresponding to its states and we do not dis- 
tinguish between sites corresponding to the same atom 
and sites corresponding to different atoms of the donor- 
bridge-acceptor system. Sometimes it is more convenient 
to treat the bridge as consisting of atomic groups such 
as, for example, the methyl groups CH2 instead of single 
atoms. Then any such atomic group can be represented 
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as a set of sites, each corresponding to a molecular or- 
bital for the group. For further simplification we assume 
that intraatomic hopping integrals are smaller that inter- 
atomic ones. This enables us to consider electron transfer 
between different sites as its tunneling between different 
atoms. 

Assume that sites " i" are associated with those atoms 
of the bridge which have the maximum coupling to the 
donor, and sites " j" are associated with atoms strongly 
interacting with the acceptor. Neglecting the intraatomic 
hopping, we can treat any site 'V which represents an 
atom at the entrance of the bridge as the origin of a chain 
along which an electron moves from donor to acceptor, 
and any site "j" as the end of the chain. Thus, we have 
a set of chains (pathways) for an electron moving along 
the bridge. In further analysis, we concentrate on the 
case when the chains weakly interact with each other, 
so we can consider them separately, as illustrated in the 
diagram below. 
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We also simulate the donor/acceptor subsystems as 
sets of semi-infinite tight-binding homogeneous chains. 
Each chain is attached to a site of donor/acceptor which 
can be effectively coupled to the bridge. As before, we 
assume that the chains do not interact. The adopted sim- 
ple model does not enable us to carry out a proper treat- 
ment of coupling of electrodes to the molecule at metal- 
molecule-metal junctions which is a nontrivial and intri- 
cated problem | 36ll4Cll4lLl42| . Here, however, we concen- 
trate on the analysis of electron tunneling from the donor 
to the acceptor through the bridge. It seems reasonable 
to conjecture that intrinsic features and characteristics 
of this process in the molecules with complex donor and 
acceptor subsystems do not strongly depend on details of 
coupling of the electrodes to the donor/acceptor due to 
a comparatively large size and complicated structure of 
these subsystems. 

In the following calculations we start from a tight- 
binding Hamiltonian for a single chain included into the 
bridge: 



H = Hn 



(4) 



whose matrix elements between states < k\ and \l > cor- 
responding to the k-th and Z-th sites of the chain (the 
states are supposed to be orthogonal) are given by: 



(Hq)u = akSki 



(Hi)ki = Vki , 



(5) 



where Ski is the Kronecker symbol, V&;= when k = I, 
and only states associated with valence electrons are con- 
sidered. The diagonal matrix elements ak are ionization 



energies of electrons at sites fc, and Vfc; = Vik includes 
both direct and exchange energy contributions for an 
electron to transfer between the fc-th and Z-th sites. As 
well as it was carried out in [37l l38| we also take into ac- 
count self-energy corrections arising due to the coupling 
of the donor (H-p) and acceptor {Hj) to the bridge. As 
a result, we arrive at the effective Hamiltonian for the 
chain, i.e., 



H, 



If 



H + H- 
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Assuming that the chosen chain starts at the i-th site and 
ends at the j-th site and generalizing the results of [3l| 
for the case when we have several donor and acceptor 
sites coupling to the bridge, we obtain nonzero matrix 
elements of H-p and Hj[ in the form 



(H-d)u = (£x>)i — E 



X> 



E 



A 2 



E 



(7) 



(8) 



Here, T> m i and Aj n are respectively coupling strengths 
between the m-th donor site or the n-th acceptor site 
and the i-th or j-th site of the bridge, and er m .„ = 

\ \d m ,n — i^J^lm.n " ^m,n} are self-energy corrections 
of the semi-infinite chains attached to the correspond- 
ing sites [s^. The parameters 6 m ,n = E — t m ,n] £m,n 
and 7 mjn are ionization energies of electrons at the cor- 
responding donor /acceptor sites, and nearest-neighbor 
hopping integrals for the chains. Summation in Eqs. 
(7) and (8) is carried out over all donor/acceptor sites 
coupled to the bridge. When E — e mj „ is smaller than 
7m, n, we have cr m n w —ij m ,n- In the opposite limit 
I E — e m ,„ |>> 7,„.„ this quantity takes on values close 
to zero. 

Due to the presence of the self-energy corrections the 
eigenvalues of the effective Hamiltonian (6) include imag- 
inary parts which represent broadening of the bridge en- 
ergy levels Ei originating from the coupling of the bridge 
to the donor and acceptor systems. The energy levels 
are further broadened if we take into account scattering 
processes in the bridge. 

Treating the long-range ET processes, we can neglect 
the direct coupling of donor sites to acceptor sites. The 
probability amplitude for the transition between the ru- 
th donor site and the n-th acceptor site is given by (2), 
and Gij is the matrix element of the Green's function 
corresponding to H e f f . An electric tunneling current / 
flowing from donor to acceptor through the bridge in the 
presence of a small applied voltage V has the form: 
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where f(E) is the Fermi function, — e is the charge of 
the electron. Assuming that the applied potential varies 
linearly in the molecule, we can use the approximation 
V = VqLi,/ L mo i, where Lf, and L mo \ are the lengths of 
the bridge and the whole molecule, respectively, and Vq 
is a voltage applied across the whole molecule. Starting 
from this expression (9), and following a usual way we 
arrive at the standard formula [8fij |: 



dE T(E)[f (E - ^) - f (E - » 2 )] 



(10) 



where the chemical potentials fj,i and fj, 2 are determined 
by the equilibrium Fermi energy of the bridge Ep and 
the effective voltage V across the bridge |3q|. i.e. 

/ii = Ef + (1 — rj) e Vj = Ef — rjeV. 

Here, the parameter 77 characterizes how the voltage V is 
divided between the two ends of the bridge. The electron 
transmission function included in the expression (10) is 
given by the formula 



T(E)=4^,A i \G ij \ 2 A j 



(11) 



h3 



The summation in (11) is carried over states < i\ at 
the entrances to the bridge and states \j > correspond- 
ing to the exits from the bridge to the acceptor subsys- 
tem, therefore contributions from all possible pathways 
are contained here. The quantities Aij are imaginary 
parts of the self-energy corrections namely: 



A, = Im(Sx))i 



Aj = Im(S^) 



.1.1 



(12) 



We see that the dependence of the electron transmis- 
sion function on energy is determined by the contri- 
butions from different donor/acceptor sites, as well as 
from the Green's function matrix elements correspond- 
ing to different chains included into the bridge subsys- 
tem. When the tunnel energy E takes on a value close to 
£ m „ or to one of the poles of the Green's function, the 
relevant term in (11) can surpass all remaining terms. 
Thus, for different values of E different donor/acceptor 
sites and different pathways can predominate in the elec- 
tron transfer through the bridge. 

This model was applied in numerical simulations of the 
electron transport in a porphyrin-nitrobenzene molecule 
013. The latter was chosen for it has complex donor 
and acceptor subsystems connected with a relatively sim- 
ple bridge. It was shown that low temperature I — V 
characteristics exhibited a step-like behavior. The steps 
originate from the structure of the electron transmission 
T(E) within the energy range where f(E— /ii) — f{E— /i 2 ) 
differs from zero. The steps disappear as the temperature 
rises. Also, the structure of I — V curves could be sig- 
nificantly affected due to electronic dephasing analyzed 
below. 



III. DEPHASING EFFECTS 

Nuclear motions existing in realistic molecules give 
rise to the electronic phase-breaking effect. This affects 
electron transmission function and the electron current 
through the molecule. When the dephasing is strong 
(e.g. within the strong thermal coupling limit), the pre- 
dominant ET mechanism is sequental hopping which re- 
places the coherent tunneling dominating at weak de- 
phasing. Typically, ET processes occur within an inter- 
mediate regime, when both coherent and incoherent con- 
tributions to the electron transmission are manifested. 
The general approach to ET studies in the presence of 
dissipation is the reduced dynamics density matrix for- 
malism (see e.g. 0, E^). This microscopic computa- 
tional approach has advantages of being capable to treat 
transition from nonadiabatic to adiabatic ET and to pro- 
vide the detailed dynamics informations. However, these 
informations are usually much redundant as far as sta- 
tionary ET processes are concerned. Therefore we choose 
an alternative approach using the scattering matrix for- 
malism and phenomenological Buttiker dephasing model 
This enables us to analytically treat the prob- 
lem, and the results agree with those obtained by means 
of more sophisticated computational methods, as it was 
demonstrated in the earlier work |27j . 

To proceed analytically, we restrict further considera- 
tion with the case of a single-site bridge. We assume that 
the bridge is coupled to a phase randomizing scatterer. 
Following the Buttikers work we simulate the scatterer 
as shown below. Here, we present the ET process by 
tunneling through two barriers separating the donor and 
acceptor from the bridge site where the dephasing could 
occur. The squares represent tunnel barriers, and the 
triangle imitates a scatterer coupling the bridge to a dis- 
sipative electron reservoir. 

This simple model could be succesfully applied to visu- 
alize and analyze various phase-breaking processes such 
as scattering on vibrionic phonon modes in the molecule 
and/or scattering on the phonons which represent the 
thermal surrounding. All information concerning char- 
acteristics of the scattering processes is contained in the 
dephasing parameter e introduced below 
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The scatterer coupled to the bridge is described with 
the 4x4 scattering matrix s relating outgoing wave 
amplitudes a' 1; a' 2l a' 3 , a' A to the incoming amplitudes 
01,02,03704- Introducing the phenomenological param- 
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eter e which measures the dephasing strength we arrive 
at the result 1291: 
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The parameter e takes on values within the interval [0, 1]. 
When e = the bridge is detached from the scatterer, and 
the ET is completely coherent. In the opposite limit e = 1 
the electron is surely scattered to the reservoir and then 
re-emitted from there, which results in complete phase 
randomization. 

Also, the electron tunneling through a single barrier 
is described with the transmission and reflection ampli- 
tudes (t and r, respectively). These are the matrix ele- 
ments of a 2 x 2 matrix: 



s e i 



t 



(14) 



Combining Eqs. (13) and (14) we arrive at the expression 
for the scattering matrix S which gives the relation of the 
amplitudes b[, b' 2 , a' 3 , a' 4 to the amplitudes b±, b%, a 3 , a^. In 
other words, this matrix describes the ET in our system 
with a sing le-site bridge. The scattering matrix S is given 

by' 



n a sing 



/ r\ + a 2 r 2 cctit-2 



S = Z~ 1 



atyti T 2 J rO?r\ a(3rit 2 f3t 2 
0t\ a(3r\t 2 /3 2 r\ ar\r 2 — a 



\ af3t\r 2 



(3t 2 



ar x r 2 



P 2 T2 ) 

(15) 

where Z = 1 — o?r\r 2 \ a. — \f\ — e; (3 = y/e; T\,i and 
ti ; 2 are the amplitude transmission and reflection coeffi- 
cients for the two barriers. The above expression could be 
generalized to include an arbitrary number of the bridge 
sites, as shown in Ref. [27j . Within this approach we 
introduce the effective electron transmission T e f f (E) : 



T, 



eff 



Ja 
Jd 



(16) 



where Jd is the incoming current from the donor and J a 
is the outgoing current from the acceptor. An electron 
could be injected into our system via four channels indi- 
cated in the above schematic drawing. Outgoing currents 
in these channels (Jj) are related to the incoming ones 
(Ji) by means of the transition matrix T (Ty — |jSy| ) : 



J'j — 2_, ' ." Ji 



(17) 



Here, Jd = Ji, J a — J 2 - To provide the charge conserva- 
tion, the net current must be zero in the channels 3 and 
4 connecting the system with the reservoir, so we have 



J3 + Ji 



(18) 



There are no grounds to anticipate any difference be- 
tween these two channels, therefore we assume J3 = J4. 
Also, we put J 2 — presuming zero income current in 
the channel 2. 

Now, we use equations (17) to express outgoing cur- 
rents in terms of the incoming ones, and we substitute 
the results into (18). This gives 



J3 — Ja — Jd 



Ki(E) 
! - R(E) 



(19) 



where 



K 1 (E)=T 31 +T 41 , 



R(E) = T 33 + T44 + T43 + T34. 



(20) 



Inserting the expressions (19) into Eq. (17) for the cur- 
rent J 2 = J a, we obtain 



T eff (E) = T(E) 



K^K^E) 
2 - R{E) 



(21) 



Here, the functions T(E), K 2 (E) and R(E) introduced 
in the expression (21) are given by: 



T(E)=T 21 , 



K 2 {E) =T 23 + T 2i . 



(22) 



To simplify further analysis we assume both barriers to 
be identical, so we have r\ = r 2 = r and ti = t 2 = t. 
Then we easily arrive at the analytical expression for the 
effective transmission function: 



Teff(E) = 



t 2 (E)(l + a 2 )[l -a 2 {l -t 2 (E))} 
2[l + a 2 (l-t 2 (E))} 2 ' 



(23) 



In the coherent limit (e = 0, a — 1) the transmission 
function T e ff(E) coincides with T(E) given by Eq. (11), 
so we can express the tunneling parameter 6(E) = t 2 (E) 
in terms of the earlier introduced Green's functions and 
self-energy corrections: 



6(E) = 



2ff 
1 + 5 



(24) 



Assuming the identical barriers, the function g(E) has 
the form 



9(E) = 



2A 



^(e-e^ + t? 



(25) 
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In this expression, Ai = A2 = A, Ai^ represent imag- 
inary parts of the self-energy corrections introduced by 
Eqs. (12), and the Green's function for the single-site 
bridge is given by 



G{E) = 



1 



E-E 1 + ?T 



(26) 



where T — 2A and E\ is the site energy. The width of this 
energy level is determined with the parameter T. Using 
these expressions, we have to keep in mind that both 
site energy E\ and the parameter T change their values 
when we attach the bridge to the dephasing reservoir. In 
our system phase breaking effect originates from nuclear 
motions which could be described as the set of phonons 
therefore E\ and T are replaced by |3fij : 

E 1 = E l + Re£ p/l (E); f = T + ImS p/l (E) (27) 

where £ p h (E) is the electronic self-energy correction due 
to the electron-phonon interaction. Nevertheless, the ex- 
pression (25) for the function g(E) maintains its form 
provided that the replacement (27) is carried out. 
So, we get 



Teff(E) 



g(E)(l + a*)[g(E)(l + ai 



I -a 2 



[g(E)(l-a* 



1 



a* 



(28) 



As it was already mentioned for the coherent tunneling 
(e = 0, a = 1) the transmission function T e ff(E) reduces 
to T{E). In the opposite limit of the totally incoherent 
transport we arrive at the expression: 



Teff(E) 



9(E) 



2A 



(E- E x f +T 2 + 2A 



(29) 



To further simplify our analysis we proceed within the 
wide-band limit. Correspondingly, we assume that the 
self-energy terms in the expression (25) do not depend on 
energy. Then the coherent electron transmission (e = 0) 
shows a sharp peak at the energy E\ (see Fig. 1) which 
gives rise to a step-like form of the volt-ampere curve 
presented in the Fig. 2. When the dephasing reservoir is 
attached to the electronic bridge (e 7^ 0) the peak in the 
transmission is eroded. The greater is the value of the 
dephasing parameter e the stronger is the erosion. When 
e takes on value 0.7 the peak in the electron transmis- 
sion functions is completely washed out as well as the 
step-like shape of the I — V curve. The latter becomes 
linear corroborating the well-known Ohmic Law for the 
sequental hopping mechanism. 

To proceed in a more realistic way we must express 
the phenomenological dephasing strength e in terms of 
some relevant energies characterizing the ET process. To 
obtain the desired expression for e we compare our re- 
sults with those recently presented in Ref. [35j. In that 
work the correction SI to the coherent tunnel current 
via a single-site bridge is calculated using noncquilibrium 
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FIG. 1: The effective transmission function (36) versus en- 
ergy. Calculation are carried out for a single electronic state 
bridge at T = 70K, E x = -0.4eV, Ai = A 2 = lOmeV. 
Curves are plotted assuming that the dephasing parameter e 
accepts values 0.0, 0.1 0.3 and 0.7 (from the top the bottom). 




0.5 1 1.5 
Voltaae (V) 

FIG. 2: Dephasing effect on the current-voltage characteris- 
tics. The curves are plotted for e = 0.0, 0.1, 0.3 and 0.7 from 
the top to the bottom. The effective electron transmission 
function shown in the figure 1 was used in calculations of the 
current. 



Green's functions formalism. The relevant result (See Eq. 
(33) in Ref. |23) is derived in the limit of weak electron- 
phonon interaction when T p h << T (T p h — ImEph). It is 
natural to assume that the dephasing strength e is small 
within this limit, so we expand our expression for T e ff(E) 
in powers of e. Keeping two first terms in this expansion, 
we get 



T eff {E)^g\E)(l+ e -[l-2g\E)] 



(30) 



We employ this approximation to calculate the current 
through the bridge, and we arrive at the following ex- 
pression for 81 : 
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SI 



zr / dE- 



cA 



r + T ph 



Pel(E) 



(E - Ei) 2 - 4A 2 
(E-E^ 2 +4A 2 



x [/(-B — A*i) — /(-E — /x 2 )] . (31) 
Here, p e i (E) is the electron density of states at the bridge: 



p e l(E) = ImG(E) 

TT 



(32) 



Comparing the expression (31) with the correspond- 
ing result of Ref. |35j we find out that these two are 
consistent, and we get 



r + r 



pit 



(33) 



When A >> T ph (e << 1) the bridge coupling to 
the dephasing reservoir is weak, and the superexchange 
mechanism of ET is predominating. The opposite limit 
A << Tph (e ~ 1) corresponds to the completely incoher- 
ent phonon assisted ET. The phonon contribution to the 
electron linewidth depends on energy E. Starting from 
the general expression for T p h (E) derived in Ref. |23j , we 
can write out the approximation appropriate at low tem- 
peratures within the relevant energy range [ii < E < fi% , 
namely: 

Tp h (E) = nM 2 lj dupp h (Lo) Pe i(E + Hlo) 



diopp h {uj)p e i{E - Tilj) \, (34) 



Here, M is the electron-phonon coupling strength, and 
Pph(w) is the phonon density of states. It includes con- 
tributions from phonons representing the thermal envi- 
ronment (stochastic nuclear motions) as well as from vi- 
brational modes. 

Intensity of stochastic nuclear motions strongly de- 
pends on temperature. Therefore we may expect vibra- 
tional modes contributions to predominate in the phonon 
density of states at low temperatures. Assuming a single 
vibrational mode with the frequency to exist in the 
system, we have: 



Pph(u) 



1 
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nh {w - Q) 2 + 7 2 



(35) 



In typical situations the linewidth 7 is much smaller than 
35], so p p h(oj) exhibits rather sharp maximum at u> = 
fl. We can roughly estimate T(E) approximating p p h(oj) 
given by Eq. (35) as 5(u> — 57). Then we obtain: 

T ph {E) nM 2 {p el (E + nn)9(p 1 -hn- E) 

+ pei (E - %Q)6{E - p 2 - hCl) } (36) 



where 9{x) is the step function. 

So we see that the peak in the frequency dependence 
of pph is reflected in the energy dependencies of T p h 
and the dephasing parameter e. Being taken into ac- 
count, this energy dependence of the dephasing strength 
may bring some features (peaks and dips) in the first 
and second derivatives of the current through the bridge 
/ with respect to the applied voltage V. Such features 
were intensely studied in recent papers (see e.g. Refs. 
HEHEQll)' an d we do not discuss the issue here. Fur- 
ther more, the dependence of the dephasing parameter 
e of energy usually does not bring significant changes 
in the shapes of / — V curves themselves, compared to 
those calculated assuming e to be a constant. Replac- 
ing p e i{E ± Ml) with their maximum values we obtain 
that T p h < M 2 /A, and the dephasing parameter e will 
take on small values provided that M < A. Therefore, 
we have grounds to expect the dephasing rate to be re- 
duced at low temperatures and the steps in the I — V 
characteristics to be distinguishable. 



IV. RESULTS AND DISCUSSION 

It follows from the obtained results that the elec- 
tron transmission function contains important informa- 
tion concerning intramolecular electron transfer. Un- 
der low temperatures and moderate electronic dephasing 
T e ff(E) can exhibit series of peaks. Their location is 
determined with contributions from donor/acceptor sites 
predominating the donor/acceptor coupling to the bridge 
at a given interval of tunneling energy E, and with the 
energy spectrum of a chain of sites connecting them. Dif- 
ferent bridge chains (pathways) contribute to the electron 
transmission function to different extent, but we can ex- 
pect the primary pathways to predominate, and to pro- 
duce the most distinguishable series of peaks in the struc- 
ture of the electron transmission function in appropriate 
energy ranges. 

The structure of T e ff(E) can be revealed in the low 
temperature current-voltage characteristics for the elec- 
tronic transport through the molecular bridges. At low 
temperatures current-voltage curves can exhibit a step- 
like behavior. This originates from a step-like character 
of Fermi distribution functions at low temperatures along 
with the "comblike" structure of the electron transmis- 
sion function T e ff(E). At a given voltage V the differ- 
ence of Fermi functions in the Eq. Hl()|) takes on nonzero 
values only in the interior of a certain energy range in- 
cluding Ep. Within a linear approximation the size of the 
range SE is proportional to the magnitude of the voltage 
V applied across the bridge. Therefore the magnitude of 
the current / at a given voltage is determined by the con- 
tributions of peaks of the electronic transmission function 
located in this energy range. When the applied voltage 
increases, this enhances the width of the energy interval 
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5E. 

However, the enhancement of 8E does not immediately 
give rise to an enhancement of the tunneling current I. 
The latter abruptly changes when an extra peak of the 
electron transmission function appears in the interior of 
the energy range where f{E— fix) — f(E— ^2) differs from 
zero. This explains a step-like behavior of the tunneling 
current as a function of the applied voltage. We also see 
that the steps originate from the structure of T e f f (E) in 
the foresaid energy interval. Widths of the plateaus are 
equal to the distances between adjacent peaks in T e ff(E) 
and magnitudes of sudden changes in the current between 
the plateaus correspond to the heights and shape of these 
peaks. 

At higher temperatures Fermi distribution functions in 
the integrand of the Eq. (|1U|) lose their step-like character 
and the plateaus are washed out. Another reason for the 
structure of T e ff(E) (and I — V curves) to be eroded is 
the electronic phase-breaking effect which arises in com- 
plex molecules due to nuclear motions. We expect, how- 
ever, that at low temperatures the electronic dephasing 
effects could be reduced so that the structure of T e f f (E) 
could be revealed. 

On these grounds we believe that experiments on elec- 
tronic transport through macromolecules at low temper- 
atures (T ~ IK) could be useful to obtain an additional 
information upon characteristics of the intramolecular 
electron transfer. Namely, comparing the structure of the 
electron transmission function reconstructed on the basis 
of experimental I —V curves with that obtained as a re- 
sult of calculations, we can make conclusions concerning 
actual primary pathways of electrons through intramolec- 
ular bridges, as well as sites of the donor/acceptor subsys- 
tems involved in the ET process at different values of the 
tunnel energy. The above comparison can be succesfully 
performed for the molecules including bridges which can 
be reasonably well simulated as sets of weakly interacting 
chains. 

The simple model of donor-bridge-acceptor system 
adopted in the present work can provide us with the 
results suitable for a quantitative comparison with the 
results of proposed experiments even avoiding a proper 
and reliable calculation of some parameters such as the 
equilibrium Fermi energy of the bridge, and the effective 
voltage V applied across the bridge. The lack of infor- 
mation about proper values of Ep and V produces an 
uncertainty in the location of the origin at the "voltage" 
axis so we cannot identify steps of the I — V curves sep- 
arately. Nevertheless, changes in the values of Ep and 
V do not influence the electron transmission function for 
a given molecule, therefore the structure of series of the 
peaks remains hxed. This enables us to identify some se- 
ries of peaks analyzing sequences of widths of successive 
steps of the current-voltage curves. Such analysis also 
can give reasonable estimations for the Ep and V for the 
chosen molecule. 



In summary, the purpose of the present work is to 
analyze how characteristics of the low-temperature elec- 
tronic transport through macromolecules can be used to 
get an additional information concerning long range elec- 
tron transfer. It is shown that at low temperatures and 
moderate electron dephasing the electron transmission 
function reveales a "comblike" series of peaks which con- 
tain information about donor/acceptor sites effectively 
participating in the ET processes and primary pathways 
of electrons tunneling through the bridge. This impor- 
tant information can be obtained as a result of analysis 
of experimental low temperature current- voltage charac- 
teristics for chosen macromolecules. 
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